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Origin of extragalactic 7-ray background (EGRB) is still a matter of debate. EGRB can either have truly 
diffuse or unresolved discrete point sources origin. Majority of the Fermi and EGRET detected identified 
sources are blazar. So, they are expected to be a significant contributors to the EGRB. In order to estimate 
their contribution to the EGRB one needs to construct their luminosity functions. Here, we examine and 
construct the evolution and luminosity function of blazars using first LAT AGN catalog. We consider both pure 
luminosity and pure density evolution models for flat spectrum radio quasars (FSRQs). We also describe the 
methodology to estimate the contribution from misaligned blazars to the EGRB. 



1. Introduction 

Observational 7-ray astronomy got a massive boost 
after the launch of the Fermi Gamma-ray Space Tele- 
scope (Fermi) on 11th June, 2008. The Large Area 
Telescope (LAT) onboard Fermi provides an increase 
in sensitivity by more than an order of magnitude over 
EGRET and AGILE. 

The origin of extragalactic 7-ray background 
(EGRB) is one of the fundamental unsolved prob- 
lem in astrophysics. EGRB can arise from some dif- 
fuse processes like black hole evaporation, large scale 
structure formation, matter- antimatter annihilation, 
etc (e.g., [H, 0, Alternatively, due to the limited 
instrument sensitivity, unresolved 7-ray sources could 
contribute significantly to the observed EGRB. 

Recently, the EGRB has been rederived from Fermi 
Large Area Telescope (LAT) data 4]. EGRB is consis- 
tent with a featureless power law of index (2.41±0.05). 
The integrated flux above 100 MeV is (1.03 ± 0.17) x 
10"^ photon cm^^ s^^ sr^^. 

Fermi already detected ~ 1500 sources in 11 months 
of observation (First Fermi catalog; [Hi) and is ex- 
pected to detect many more. In contrast to the less 
than 100 active galactic nuclei (AGNs) in the EGRET 
catalog, the first Fermi-LAT AGN catalog @ con- 
sist of ~ 600 blazars. Fermi also detected 11 non- 
blazar AGN 0, three normal (M 31, LMC and Milky 
Way) and two starburst galaxies (M 82 and NGC 253). 
Thus, Fermi provides an excellent opportunity to un- 
derstand the origin of EGRB. 

Since the dominant class of identified EGRET and 
Fermi sources are blazars, they are expected to con- 



tribute to the EGRB. The contribution from blazars 
to the EGRB is estimated using two different ap- 
proaches: 

(a) Considering a relationship between the blazar 7- 
ray luminosity with their luminosity at some other 
wavelength, their 7-ray luminosity function is as- 
sumed to be the scaled luminosity function at that 
wavelength (e.g. , [1, i [IS [n!, [H, [H, [H, [H [l|] ) . 

(b) Constructing the 7-ray luminosity function di- 
rectly from observed 7-ray blazars (e.g., [l3, [H, [H, 

In an earlier work [13, [3 following a similar ap- 
proach of Chiang & Mukherjee [l^ but with almost 
twice the number of sources, we constructed the 7-ray 
luminosity functions of flat spectrum radio quasars 
(FSRQs) and BL Lacs from EGRET-detected sources 
only. Wc found strong positive evolution for FSRQs, 
while BL Lacs did not show any significant evolu- 
tion. The maximum contribution from blazars to 
EGRB was found to be ~ 20%. However, due to 
the limited number of sources in the EGRET catalog, 
the luminosity function and hence, the contribution 
from blazars was not well constrained. Miicke & Pohl 
[2H estimated the contribution from FSRQs and BL 
Lacs to the EGRB in the context of AGN unification 
paradigm, and found that unresolved blazars (both 
"ahgned" and "misaligned") contribute 20% - 40% 
to the EGRB. In a similar approach, Dermer 
estimated that FSRQs and BL Lacs contribute ~ 
10% - 15% and 2% - 3% at 1 GeV respectively. From 
the first three months of Fermi observation, Abdo 
et al. [2^ found that FSRQs exhibit strong evolu- 
tion while BL Lacs show no evolution. Using the 1st 
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Fermi catalog, Abdo et al. [2J] found from the source 
count distribution that the point source contribution 
is ^ 16% of the EGRB. 

Here, we construct the luminosity function of FS- 
RQs and BL Lacs following the methodology de- 
scribed in Bhattacharya, Sreekumar & Mukherjee 
[l7j . Utilizing the Fermi observed sources we con- 
struct a complete source list, apply the V/Vmax 
method to search for the presence of evolution, find 
parameters by modified V/Vmax method, and fi- 
nally construct the luminosity function using 1/Vmax 
method and maximum likelihood estimator. 

In the next section we discuss the evolution of blazar 
source class. We discuss in section 3 the luminosity 
function construction of FSRQs. The methodology to 
calculate the misaligned blazars contribution to the 
EGRB is given in section 4. We conclude with our 
discussion in section 5. 



2. Evolution of blazars 

We study the luminosity evolution of blazars by 
considering a sample that consists of sources from first 
LAT AGN catalog having ha or above detection sig- 
nificance with 7-ray flux (E > 100 MeV) > 5 x 10~* 
ph cm~^s~^. Abdo et al. f2l| showed that Fermi- 
LAT detects spectrally hard sources at flux level gen- 
erally fainter than those of soft sources. They showed 
that the LAT sample (including all sources) is com- 
plete above 7 x 10~^ ph cm~^s~^ (Figure 1 of Abdo 
et al. (24I). However, the maximum photon index 
for blazars is 2.98, much lower than the maximum 
photon index of all sources of 1st Fermi catalog, and 
from figure 1 of Abdo et al. [2^ , one can see that the 
blazar sample is complete to a flux below 7 x 10~* 
ph cm^^s^^. As a secondary check to the complete- 
ness of our sample above 5 x 10~® ph cm~^s~^, we 
performed the < V/Vmax > test of this sample and 
obtained < V/Vmax > that matches within la of a 
sample of limiting flux 7x 10~* ph cm~^s~^. Our sam- 
ple consists of 118 FSRQs. and 39 BL Lacs (25 with 
known z). The < V/Vmax >FSRQ comes out to be 
0.63 ± 0.03, which indicates strong evolution. For BL 
Lacs (with known z sources only), < V/Vmax >BLLac 
comes out to be 0.49 ± 0.06, which indicates no sig- 
nificant evolution. If we include the sources with 
unknown z considering they are at the average z of 
the sample, < V/Vmax >BLLac comes out to be 
0.55 ± 0.05. Since, a good fraction of BL Lacs in our 
sample do not have redshift information, it is not pos- 
sible to construct a meaningful luminosity function of 
BL Lacs at this moment. 

We consider both pure luminosity evolution (PLE) 
and pure density evolution (PDF) for FSRQs. 

Both power law ((1 -I- z)^) and exponential 
(exp(T(z)/r)) evolution function have been consid- 
ered. Here T{z) denotes look-back time. The mod- 
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Figure 1: Upper panel: The variation of < 



> with 



evolution parameter /3. Lower Panel: The variation of KS 
test probability with /3. 



ified < V/Vmax > method is used to find the evolu- 
tion parameter. For the optimum parameter value, 
V/Vmax is expected to be uniformly distributed be- 
tween and 1, thus < V/Vmax >= 0.5. We find 
TpLE = 0-2lto^o2- For each value of evolution param- 
eter (r), the distribution of V/Vmax is compared with 
a uniform distribution using KS-test which shows dis- 



0.05 



tribution of V/Vmax is uniform for tple — 0-21 _o 
The optimum value of Pple is found to be 2.11q4. 
Similarly, for pure density evolution, the parameter 
values of exponential and power law evolution models 
are, tpde = Q.Uto^^l and f3pDE = S.I^q^s respec- 
tively. The upper panel of Fig [1] shows the variation 
of < V/Vmax > with different values of power law 
PDF parameter fipoE- Horizontal dotted lines show 
the la error in < V/Vmax >■ 

We derived the average 7-ray photon index (s-y) fol- 
lowing the methodology given in [13, [IB] • The values 
of < > for FSRQs and BL Lacs are 2.45 ±0.16 and 
2.17 ± 0.22 respectively. 



3. Luminosity function construction of 
FSRQs 

For pure luminosity evolution, after de-evolution of 
the sources at z = 0, we constructed their luminos- 
ity function by 1/Vmax method which suggests a lu- 
minosity function described by a broken power law. 
We then used maximum Likelihood function for the 
redshift distribution of sources, considering a broken 
power law luminosity function. 



Lq < Lb ■ 



/^o X I ^ j Lo> Lb ■ 



(1) 
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Figure 2: Likelihood parameters contour plot of 
luminosity function for power law PDE model. Upper 
Panel: 68.3%, 95.4% and 99% {MogC = 1.15,3.08,4.61 
respectively) likelihood contour plot (in to out) for ai. 
Lb (left), cti (middle), a2 and a2 (right), Lb. Lower 
Panel: variation of negative log likelihood values for Lb 
(left), ai (middle) and a2 (right), the intersections of the 
solid line with the curve represents the la value. 



Figure 3: Distribution of FSRQs with redshift. 
Histogram shows the Fermi observed distribution. The 
red curves represent PLE models. The black curves 
represent PDE models, solid and dash curves are 
exponential and power law models respectively. 

4. Misaligned blazar contribution 



Table I Luminosity Function Parameters 





Lb (in lO"*' erg s^^) 


Ql 


Q2 


TPLE 




1.5 ±0.3 


*-O-0.3 


Pple 




1.5 ±0.2 


7 ,+2.4 
'•1-1.1 


TPDE 


61 03+21-68 


2.0 ±0.1 


7+0.4 
' -0.2 


Ppde 


on. 1 0_g Q2 


1.8 ±0.1 


^■'-'-0.2 



Here 4>{Lq) is the de-evolved luminosity function and 
00 is the normalization of the luminosity function. 
The break luminosity (Lb) and the lower and upper 
end index (ai and a2 respectively) are constrained by 
maximizing the likelihood function. 

For pure density evolution, we consider a similar 
form of luminosity function and constrained the pa- 
rameters by maximizing the Likelihood function for 
the redshift distribution of sources. 

The derived luminosity function parameters values 
are given in Table HI Fig [2] shows the contour plots 
of different luminosity function parameters for power 
law PDE model. 

The normalization constant of luminosity function 
has been determined by equating the total number 
of FSRQs in our sample to that expected from the 
model luminosity function. The source distribution 
(dN/dz) has been over plotted with the model distri- 
bution for PLE and PDE models in Fig [31 Exponen- 
tial and power law evolution models are considered for 
both PLE and PDE models. 



According to the AGN Unification scenario, FR II 
galaxies are the parent population of FSRQs whereas, 
FR I galaxies are the parent population of BL Lacs. 
The AGN jet emission is expected to fall very rapidly 
with increasing jet to line-of-sight angle. Nevertheless, 
considering the large population of these misaligned 
blazars compared to nearly-aligned ones, these sources 
could contribute significantly to EGRB. Fermi has al- 
ready detected 11 misaligned blazars. 

The observed jet luminosity L{d) at a jet to line-of- 
sight angle {9) can be written as 

m - m X m (2) 

From the knowledge of £^{9) one can in principle, es- 
timate the contribution from off-axis AGNs to the 
EGRB. 

4.1. Construction of ^{9) 

Our aim is to model the SED of all Fermi-detected 
off-axis sources using standard leptonic models with 
structured jets. From the knowledge of these SED 
modeling and VLBI observations one can fix the emis- 
sion processes in the jet and other jet parameters and 
can construct S,{9). Next step is to constrain ^(0) 
by comparing 7-ray detected misaligned blazars with 
equivalent on-axis objects (by matching sources with 
similar jet power, F) and then validate our predictions 
by modeling few misaligned blazars that are not de- 
tected by Fermi. 

From the observed 7-ray luminosities of off-axis 
sources and their jet inclination angle one can also 
construct £,{9). The difference in the observed 7-ray 
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Figure 4: SED modeling of misaligned blazars 3C 120 



showed that the redshift of the peak in the number 
density of FSRQs is luminosity dependent. Hence, 
they considered himinosity dependent density evolu- 
tion for FSRQs (similar to that of radio-quiet AGN). 
They considered a much lower limiting flux (~ 10~^ 
ph cm~^ ^~^) s-nd hence, a much larger sample (186 
FSRQs) than us. 

For misaligned blazars, we plan for SED fitting of 
all detected off- axis sources and then construct ^(0) 
from model. This work is under progress and will be 
reported elsewhere. 
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luminosities between two radio galaxies can arise from 
two different reasons: 

a) due to their different jet to line-of-sight angle, 

b) due to their different intrinsic luminosities. 

If one considers that the total radio luminosity is 
a tracer of unbeamed luminosity, then after an ap- 
propriate scaling, one can construct ^(0) from scaled 
7-ray luminosities and jet inclination angles of off-axis 
sources. 

We have initiated this work. Our aim is to analyze 
Fermi data to date on radio galaxies and model it to 
derive jet parameters. We have completed the analysis 
and the SED modeling of 3C111 and 3C120 (Figureg] 
shows the SED fitting of 3C 120.) 

Both the sources are well fitted with simple SSC 
emission model. F is also ~5 for both the sources. 



5. Discussion 

Using first Fermi catalog, we investigated the lumi- 
nosity function and nature of evolution of FSRQs and 
BL Lacs separately in 7-rays. The nature of evolution 
of FSRQs and BL Lacs as found in this work, are in 
agreement with earlier works (e.g., [l^iliPl). We con- 
sidered both PLE and PDF models for FSRQs. Both 
power law and exponential evolution functions have 
been considered. The luminosity function is consid- 
ered to be of a broken power law form and the param- 
eters are found using likelihood analysis. From this 
analysis it is not possible to strongly conclude on the 
true evolution model. Fig |3] suggests that PDF model 
explains observations marginally better than that of 
PLE model. In contrast to earlier works of deter- 
mination of luminosity function from 7-ray observa- 
tions [13, our derived luminosity functioir^s 
well constrained. In a recent work, Ajello et al. i26| 
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